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O ■ Abstract 
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O \ By using the background field method of QCD in a path integral approach, 

we derive the equation of motion for the classical chromofield and for the 
On ' gluon in a system containing the gluon and the classical chromofield simul- 

taneously. This inhomogeneous field equation contains a current term, which 
CN . is the expectation value of a composite operator including linear, square and 

^*~i \ cubic terms of the gluon field. We also derive identities which the current 

t~^ \ should obey from the gauge invariance. We calculate the current at the lead- 

^^ ' ing order where the current induced by the gluon is opposite in sign to that 

(^ ■ induced by the quark. This is just the feature of the non-Abelian gauge field 

^P . theory which has asymptotic freedom. Physically, the induced current can be 

treated as the 'displacement' current in the polarized vacuum, and its effect 
is equivalent to redefining the field and the coupling constant. 

(D . PACS: 12.38.-t,12.38.Aw,11.15.-q,12.38.Mh 



I. INTRODUCTION 

The behavior of electric charged particles in a classical electromagnetic field has been 
extensively studied in the literature. However, the behavior of color charged particles (quarks 
and gluons) in a classical chromofield has been less studied. The gluon case is particularly 
interesting because there is no analogue in the electromagnetic theory, since photons do 
not interact with the classical electromagnetic field. One important issue is to derive the 
equation of motion for the gluon in the presence of chromofield and to define and calculate 
the classical current produced by the gluon. As we know it is easy to define and understand 
the color current for quarks in the classical chromofield. But strictly defining the color 
current for gluons classically is not an easy task, mainly due to its quantum and bosonic 
nature. We plan to address some aspects of these issues in this paper. 

At first sight, this problem looks as if it were only of theoretical interest. However, 
it turns out to be applicable to a non-equilibrium system of quarks and gluons, which is 
expected to be produced in ultra-relativistic heavy ion collisions in RHIC and LHC where a 
classical chromofield may be produced. The dynamics of a quantum non-equilibrium many- 
body system is usually studied in the closed-time-path formulation with non-local source 
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terms [|T|. This formalism needs to be extended by incorporating classical cliromofield. 
The evolution of partons can be obtained by solving the Schwinger-Dyson equation defined 
on a closed-time-path. This is formidable task because it involves a non-local non-linear 
integrodifferential equation and because of its quantum nature. There are two scales in 
the system: the quantum (microscopic) scale and the statistical-kinetic (macroscopic) one. 
When the statistical-kinetic scale is much larger than the quantum one, the Schwinger-Dyson 
equation may be recasted into a much simpler form of the kinetic Boltzmann equation by a 
gradient expansion. The Boltzmann equation describes evolution of the particle distribution 
function in the phase space of momentum-coordinate, and can be solved numerically for 
practical purpose. For a non-equilibrium system of quarks and gluons in a chromofield, 
the distribution function also depends on the classical color charge of the parton, since the 
color is exchanged between the chromofield and partons and among partons themselves. In 
this case, the Boltzmann equation also describes the evolution of the parton distribution 
function in the color space 0,^- While the Boltzmann equation describes the kinetic and 
color evolution of the hard parton system, soft partons are normally treated as a coherent 
classical field whose evolution may be described by a equation which is similar to Yang- 
Mills equation. Therefore, one should study the transport problem for hard partons in the 
presence of a classical background field. For example, in high energy heavy-ion collisions, 
minijets (which are hard partons) are initially produced and then propagate in a classical 
chromofield created by the soft partons ^^- In this situation it is necessary to derive the 
equation of the gluon and the classical chromofield to study the formation and equilibration 
of quark-gluon plasma. Hence we see that the issues we want to address in this paper have 
practical implication in ultra-relativistic heavy ion collisions. 

Since conventional QCD does not contain classical chromofield, to include the classical 
chromofield in a proper way, we follow the background field method of QCD (BG-QCD). 
This method was first introduced by DeWitt and 't Hooft 0J^. The advantage of BG-QCD 
is that it is formulated in an explicit gauge invariance manner (see below). This method 
is widely used to study confinement, gravity and hot quark-gluon-plasma, etc.. In this 
paper, we derive the equation of motion for the gluon and the classical background field to 
define the induced current in a path integral approach which fully incorporate non-Abelian 
gauge symmetry. Also some properties of the induced current are given in the framework of 
BG-QCD. 

The paper is organized as follows. In section II, we briefly introduce basics of BG-QCD 
and the Ward identity which we will use later. In section III, equations of motion for the 
classical and the gluon field are derived, and the current induced by the gluon is defined. We 
study its property and derive the Ward identity that the current obeys in the next section. 
Section V presents the result for the current at the leading order. The physical implication 
for the result is also discussed. The summary and conclusion is presented in the last section. 

II. BACKGROUND FIELD METHOD 

We know that any physical quantity calculated in the conventional QCD is gauge invari- 
ant and independent of the particular gauge chosen. However, there is still a problem of 
gauge invariance which is evident in the off shell Green's functions and counter terms. This 
problem arises because one must fix a gauge and introduce the gauge-fixing and ghost term 



to quantize the theory. The classical Lagrangian for the gauge field has explicit gauge is 
invariant while the gauge-fixing and ghost term is not gauge invariant although it is invari- 
ant under the BRS transformation. The background field method is such a technique which 
allows one to fix a gauge (background-field gauge) without losing explicit gauge invariance 
which is present in the classical level of the gauge field theory. 

By writing the conventional gluon field into the sum of a classical background part and 
the quantum one, the action in the background field gauge is given by [0-0: 



O — Oq + ^fix + J ghost + ^ sources 

Sf^x = -i^jdy{D%A]Q^^^)\ (1) 

Sghost = I d^xC'D%A]D^^f[A + Q]C\ 

Ssource = I d^x{j;Q>^^^ + fC' + C^e)- 

The generating functional for Green's function is 

Z[A, J,e,e] = J[dQ][dC][dC] exp{iS). (2) 

The variables are defined as follows: A^ is the classical background field; Q* is the gluon 
field, the variable of integration in the functional integral; C*, C are ghost and anti-ghost 
fields; ^ , ^* are external sources coupling to ghost and anti-ghost respectively. -D^-'[^] is the 
covariant derivative in the adjoint representation whose generators are {f^Y^ = if^"'-' and it 
is defined as D'^[A\ = d^5'^ + gf''^ Al= 3^5'^ -ig{t\y^ A^- F^^[A] is the field strength tensor 
for the background field and given by -F^y[yl] = d^A^^ — dyA^^ -\- gf^^'^A^^A'^. 

The gauge transformation for all fields and sources which leaves Z invariant is given by: 

A'^ = UA,U~' + IUd,U-\ 

Q'^ = UQ,U~\_J'^ = UJ,U-\ 

C; = UCU-\ e' = UiU'\ ^ ^ 

c' = ucu-\ i' = uiu-\ 



This gauge transformation is known as Type-I transformation |Tl| . Note that the background 
field transforms in the conventional way while the gluon field, the ghost field and all external 
sources transform like a matter field. The infinitesimal changes for these fields and sources 
are: 



k 

(4) 



Here we define the gauge transformation as U{x) = exp{iguj"'{x)t'\). Roman alphabets de- 
note color indices and Greek ones denote the Lorentz indices. All fields and sources are either 
in the vector form, e.g. A^ = ^^t^, or with their color indices explicitly shown, e.g. ^4° 
It is obvious that the generating functional Z is invariant under the above transformation. 



S^, 


= D'^[A]u\ 




SQ', 


= gP'Qiu;', 


5Jl = gP'^u; 


6C'-- 


= gP'^&uj\ 


sT = gP^i'^\ 


sc" -- 


= gr^C'u\ 


se = gP'^'u;''- 



But the less obvious is the invariance of Sfix- To prove it, one should remember the fact 
that t\ is a pure imaginary Hermitian matrix and that (t^)^ = —t% therefore f/"^ = U~^. 
The generating functional for the connected Green function is given by: 

W[A,J,i,^ = -i\nZ[A,J,i,l] (5) 

Field averages can be obtained by taking derivative of W with respect to their corresponding 
sources: 

5W 5W —i 5W 

«y-w^- <^'>-^- <^>--.F (^' 

The one-particle irreducible generating functional T[A, {Q), {C), {C)] is defined by Legendre 
transformation for W: 

T[A, {Q), {C), {C)] = W[A, J,e,e] - jd'x{Jl{Q^'')+t{C^) + (^)e) (7) 

where we see that j— = |^. Since T[A, (Q), {C), (C)] is a gauge invariant functional of A, 

{Q), (C), and (C), hence T[A, (Q), {C), (C)] should be invariant under the transformation 
(1). Then we have the following Ward identity |1TT| : 

In the same way, we have the parallel Ward identity in terms of W: 

z,.i^]_+,rv^^+e^+^'^)^o. (9) 

Note that Eq.(||) and (P) are identical because -^ = |^. 

III. EQUATION OF MOTION FOR GLUON AND CLASSICAL FIELD AND 
DEFINITION OF INDUCED CURRENT 

By shifting the quantum field Q* (x), a variable of integration in the functional integral, 
by a small amount v^^{x): Q^(a;) = Q^^{x) + f^(a;), we can obtain the equation of motion or 
the Lagrange equation for the quantum field Q* . The change of the functional for connected 
Green's functions W is: 

5W = -i5\nZ = Z-^ f[dQ][dC][dC] [ d^x^vi{x)exp{iS). (10) 

The fact that W should be independent of the shifting value of the integration variable leads 
to the equation of motion for Q^. 



) = Z-' J[dQ][dC][dC] r \o^^. ^^''''^ ^ J''\x)]eM^S) = 0. (11) 



'SQU^y J'^'' '' '\ SQlix 



In the absence of the external source J* the equation of motion becomes: 

6 So SSfix , , SSghost 



where we have 



^ = ^^^'^^[^]^^'^'1^]Q^, (13) 

Note that the partial derivative operation in /} ^' '^ is toward C . We know that F^^'^[A + Q\ 
transforms like a matter field and D'^J [A + Q] is commutable with the transformation matrix. 
So the first equation ^Jf = D^J[A + Q]F'^^^'^[A + Q] transforms like a matter field in adjoint 
representation. Similarly we can easily verify that the second equation transforms like a 

CO 

matter field. But it is less apparent to verify that the ghost part ^^r* also transforms like 

a matter field. Let us verify it directly. The variation of c^°"' under the infinitesimal gauge 
transformation Eq.(|) is: 

5(%^) = -gf''6(C' b ''''''[A])C' - gf^'C' D '''''\a\6C' 

' r A^( „ flder^d, ,e^ 

(14) 



= -gt^'^igt^'^C^ b ^'^\A]uj'')C^ - gf^'^C^ D ^''''^[A]{gf'^''C^uj') 

= _g2(^f aid fake ^ fkil flde-^QJ 'jj ^"^^^ yA]C'^Uj'' 

= gf'[-gf^'^C' d'"'^'[A]C'^]uj' 

Therefore we see that g^°°' transforms like a matter field in adjoint representation. 
We can rewrite Eq.(|l^) in another form: 

Dl^[A]F-'^'^[A] = {jn (15) 

where the current jo induced by the quantum field is given by: 

j^'' = -gf^b^Di''[A]{Q'''^Q^'''') - gf'^''Q'^^{D'''''^[A]Q'''^ - Df''''^[A]Q'''^) 

_g2jidajabcQdQu,bQf.,c _ D'J [A]{D'''^''[A]Q^''^ - D'^'^'^IAIQ"''') - gp^QlF'^i'^^lA] 

(16) 

The last two terms in jo are just — (5(S'/j^ + Sghost)/5Q- Obviously jo transforms like a matter 
field. 

It is easy to check 6S0/6Q = 6S0/6A. But it takes some steps to prove: 

/(S/.. + V..)),( ^(WV»..) )^ (17) 



,501 ' ' M 



M 



We start with Z[A, 0, 0, 0] where in Eq.(^ the sources are set to zero. We shift A by 6A and 
Q by —6A so that A + Q does not change. Sq is invariant in this operation. The variation 
of Z[y4, 0, 0, 0] is caused by that of (S/j^ + S ghost)'- 



5Z[A,0,0,0] = i J[dQ][dC][dC]{^-^^^^^^^^i^ - ^i^£iitf£^^}5^i exp(i^) 



= iJ[dQ][dC][dC]{-^D'^^'^[A + Q]D'^^^''[A]Qt (18) 

where S = Sq + Sfix + Sghost- We expect that it would leave Z[y4, 0,0,0] unchanged i.e. 
5Z[y4, 0, 0, 0]/(5y4^ = 0, since nothing happens for So[A + Q]. But it is non-trivial to prove 
this. We provide here a proof based on Ref. |jll| and |1^. Consider a kind of special 
infinitesimal gauge transformation defined by: SQ"^ = D^^[A + Q]u!^ where the infinitesimal 
gauge group parameter uj^{x) is defined by: 

u^{x) = M-''^\x,y)D'J[A + Q]{y)SAi{y) (19) 

where the repetition of indices implies the summation or the integration; The matrix M is 
just the one which couples to the ghost field: 

M^^{x,y)=Df[A]D'j[A + Q]{x)6{x-y) (20) 

we have following properties: 

D^^[A]{x)6Q^^{x) = M^\x,y)uj\y) = D^\A + Q]{x)6Aj^{x) (21) 

We say this infinitesimal gauge transformation is special because it can make gauge condition 
invariant in some circumstances. In this paper, we imply the covariant gauge constraint 
6{D^ri[A]Q'^'^ — a*) in the generating functional Z. If we increase a* by a small amount 6a^ 
and simultaneously make the above gauge transformation where the small parameter uj^{y) 
is given the value so that M'^^{x,y)ui^{y) = 6a^, the gauge constraint remains the same. In 
fact the new gauge constraint 

6iD^^[A]Q'^'^ - a'O = SiD"^[A]Q^'^ + D^M^^Qj, - a^ - 5a^) 

= 6{D'^[A]Q''^^ - a') ^ > 

is unchanged. The other part of the generating functional Z can be proved to be invariant 
under the above gauge transformation constrained by M'^^{x^y){J^{y) = 6a^. Hence Z is 
independent of a*. One can therefore multiplies an arbitrary function of a* and integrate 
over a* without changing the result up to some irrelevant normalization constant. Usually 
this function of a* is chosen to be a Gaussian function exp(— 'i/(i^x(a*(x))^/2a) which gives 
the gauge fixing term Sfix- One can prove that the variation 6Z caused by the above defined 
gauge transformation constrained by M'^^{x,y)ui^{y) = Sa^ is just 5Z[A, 0,0,0] in Eq.(|18D 
which is zero. Therefore we obtain {S{Sfix + Sghost)/^ A) = {6{Sfix + Sghost) /^Q)- Together 
with 5Sq/5Q = 6S0/6A, we have 

{6{So + Sf,x + Sghost)/ 5A) = {5{S^ + Sf,x + Sghost)/ 5Q) = (23) 

From {5{Sq + Sfix + Sghost)/ 5 A) = 0, we obtain: 

D'J[A]F''^''^[A] = (j^'^), (24) 

where j is defined as: 



-gf'^C^ D ''''"^[A]C^ - gp^C^D>'^''\A + Q]C^ + lgf'i^Q>'^''Dl\A]Q''^^ 

(25) 

where the last three terms are just S{Sfix + S ghost) /^A'-^. Because of the identity 
D'^[A]Dl''[A]F''f''''[A] = 0, we have D'^[A]{j^''^) = 0. 

Eq. (|T5|) , Eq.([T7|) and Eq. (|2^) are our main results. They guarantee the consistent def- 
inition of the induced current. As we see it is easier and clearer to prove the equation of 
motion for (Q), (|^) = 0. But it is less apparent for (|^) = 0, the equation of motion for the 
classical background field, in a quantum field formalism. As a contrast, it is straightforward 
to obtain the Lagrange equation of motion for a field in the classical field theory, but the 
current situation is subtle because the treatment of the quantum field, the gauge fixing term 
as well as the ghost field and is highly non-trivial. The proof for the equivalence of two 
equations is guaranteed by the intrinsic relation between the the gauge fixing term and the 
ghost sector, or in other words, the equivalence and the consistency of two equations is a 
natural outcome of the symmetry behind the non-Abelian gauge field theory. We notice 
that Ref. [|1^ also derives Eq.(p4D from a canonical formalism, but what we use here is a 



path integral approach where the symmetry of non-Abelian gauge field is explicitly fulfilled 
and exploited. 

IV. IDENTITY FOR INDUCED CURRENT 

We now try to derive the equation for the current from the Ward identity. In order to 
derive the equation for the current from the Ward identity, we add to Ssource a new source 
term / d'^xY^j'^''^ where the current j induced by the quantum field couples to its external 
source Y. Since j transforms like a matter field, we require that Y also transform in the 
same way to ensure the gauge invariance of the new source term. 

Before we proceed, some points in the derivation should be clarified. We distinguish two 
cases: 

Case I. There is no constraint on values of (Qp = SW/6J'^'\ (C*) = SW/d^"- and 

(C ) = —6W/6^^ in intermediate steps. Therefore, during the derivation there is no relation 
between external sources (J, ^, ^) and the background field A. Hence sources and the 
background field are independent variables. Because we do not impose any specific value on 
(j) in the whole process, the source Y remains irrelevant of A. 

Case II. (Q), (C) and (C) are set to zero in each step, as assumed by many other 
authors PJTO[], hence specific relations between sources (J, ^, ^) and A are built up by 
^W[J,^,^]/SJ^''' = 0, 5W/S'^ = and 5W/5C = 0. Each source is a functional of A and 
vice versa. Through A, there are functional mappings between any two sources too. These 
sources are no longer independent of each other. But we do not assume any specific value 
for (j), so the source Y is still independent of A. In this case, we have two independent 
variables, Y and A. 

So in both cases Y is always an independent variable. We will see in either case, i.e. 
whether {Q), (C), (C) = or not, the equation for (j) is the same, though the Ward iden- 



tity for each case is different. Here we assume case I(For case II, see Appendix A). After 
introducing a new source term, the Ward identity (||) is extended to: 

5W ■•, -dW -dW -dW -i5W 

D-M\j^^ + sr'p-^ + Jijj, + e^ + e-^) - 0, (26) 

where 6W/6Ajj^ is given by: 

Substituting Eq.(|7D into Eq.(|2|) and using the identity D'^[A]Di^[A]F''^'^''[A] = 0, Eq.® 
becomes: 

, Sj"''' ... .SW -dW -SW -i5W 

- D-mr^) + D-[A]Y^'{j^^) + 9r\Y^^^ + -^^^ + ^V + ^V^ ^ °' ^''^ 

Because we want to go further, we must keep all sources in Eq.(|^). If we set all sources 
vanish at this stage, we get Dj^^[A]{j'^''') = which we obtain in the former section. Tak- 
ing derivative with respect to Y^{y) for Eq.(|^), setting all sources vanish, and using 
Dj,^[A]{j^''') = 0, we obtain: 

D^^[A]{x){T[j^^^{x)r'{y)])+z6\x-y)gf^A^^{x){^{x)) 

+igS\x - y)f'\f'\x)) + zd,[6\x - y){^{x))] = 0. ^ ^ 

Integrating Eq. (|29|) over x, and noting that the complete differential vanishes, we obtain: 

f^ I d'xA[ix){T[f^^ix)r\y)]) + ^f^A'^iy){^J^iy))+^f^ = 0. (30) 

In addition, integrating Djj'[A]{j'^'''{x)) = over x, we obtain: 

d^xf'Al{x){j^''''{x)) = (31) 



Eq. ( p9l - pl| ) are identities for j. In above formula, we use the language of functional path 
integral and the expectation value is taken in a functional sense. We can translate path 
integral expressions into the language of quantum operator due to the exact correspon- 
dence between two approaches. Then the expectation operation is carried on the inter- 
acting vacuum state, i.e. that of the full Hamiltonian H = Hq + Hi where H^ denotes 
the interacting Hamiltonian. In above formula, all operators are in the Heisenberg picture. 
According to Gell-Mann-Low's theorem, we can rewrite them in the interacting picture: 
((■ ■ ■)h) = {n\{- ■ ■)H\n) = (0|(- ■ ■)/e-*/'^*^'^W|0) where \n) is the vacuum state for the full 
Hamiltonian H and |0) is that for the free Hq, and 'H' and T denote the Heisenberg and 
interacting picture respectively. 

Note that -g^{x) and j{x) are composite operators, so they should be understood as 
normal product in order to be well defined in taking their expectation values on the vacuum 
state. If we assume, \Q) ~ |0), Eq. (p9D is then: 
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D'^[A]ix){0\T[j'^^^ix)f^'{ym = 0. (32) 

where we drop in j linear terms of Q; Also we have used (0|-^73-(x)|0) = and (0|j''''(x)|0) = 

because -fjri^) ^^^ J'^'K^) ^^^ treated as normal products. 

Because we notice that ^ transforms in a good manner, similarly we can also add to the 
action a new source term, where ^ couples to its external source, to obtain a new identity. 
For simplicity, we drop linear terms in j which results in that ^ is independent of A. Defining 

d^y^"'' = fj^, we can check the transformation rule: {!£)''"'''''' = {U (^)^V"l)'^^ We 

add to the action a new source term y°" ( Za j • '^° make the source term invariant under 
the gauge transformation, we assume that Y transforms in the same way as (^) does, i.e. 
^ui^ = {UYn,yU~^)^°' . The infinitesimal variation of Y is then: 

SY^^ = -guj\f^^Y^/ - r'^Y^:) (33) 

The Ward identity becomes: 



M L J 5Af ' ^ ^•^ 1^" ■' /^^ ^ SYb^ 



+^f^'('^^l|+e^fP+e'f) 



(34) 



Because ^ is independent of A, only |^ finally contributes to ^. Repeating the same 
procedure as in the derivation of Eq.(^) and (|30|) , we can obtain: 



D'M]i^)(nj''(^)j^(y)]) 



5A^ y-n ' J \ 5A^ 



+ig6\x - y)(f ^'(f^(x)) + /'=-(f^(x))) = 0. 
The integration form is: 



(35) 



X^f3,l XAf3,k i:-f3,a 



gf^ / rf^x4(x)(T[j^'^(x)^(y)]) + ig{r{^^{x)) + r^j^^iv))) = 0. (36) 



We can take derivative with respect to Y for Eg . (|28|) arbitrary times and get a series of 
identities, each of them shows a relation among correlation functions for j of the lower order 
and those of the higher order. This is similar to the BBKGY hierarchy in statistical physics 



V. CURRENT AT LEADING ORDER 

In this section, as an example, we calculate the expectation value of the current (j'^'*(a;)) 
at the leading order. Since j^'*(x) is a composite operator, in order to be well defined, it 
is treated as normal product. If the expectation value j'^'*(x) is taken between the vacuum 
state of a free field Lagrange, we simply have (0|j'^'*(x)|0) = 0. So {j^'^{x)) starts from the 
loop level. Here we consider the lowest order result at the one-loop level. At this level, 



the relevant term in the action is — j d'^xj^A^''^ . The corresponding term in Hi is then 
J d^xf^Af^'\ so —i J dtHi = —i J d^xf^A^'\ Hence we obtain: 

{r\^)) = -^ J d'yAi{ymT[f^^{x)fJ{y)]\0) (37) 

where j^o is part of j in Eq. ( pSj) , which is independent of A and has only square terms of 
the Q field; — / d^yAl{y)j'^o{y) is part of the classical action which contributes to the lowest 
order, j^o is given by: 






ahcr^" o ^/^c i^ „ fabcr^'' aiir^r. ^ / 



where we have not included terms proportional to Q, QQQ and QCC . Because there are 
only square and cubic terms of Q in j, linear terms of Q in the action do not contribute to the 
result. We drop the QQQ term both in j and in the classical action because it contributes 
to the two-loop level. 

Note that Eq.(|3) should satisfy Eq.(|T]): 

d'xf'^AUx){3^'\x)) = -tf f d'xd'yA'(x)Aiiy){0\T[fJ{x)fJiy)]\0) = (39) 



Actually, for the leading-order j given in Eq.(^), Eq.(pO|) and Eq.(^) are equivalent. Notic- 
ing in Eq.(|3T|) the term if^''-' A'-^{y){^^{y)) gives an additional strong coupling constant g 
relative to rest terms, we can neglect it for the leading-order j and obtain: 

f^ ld'xA^(x){T[3'J{x)3':i{y)])+ir{3^\y)) = (40) 



where only j^o contributes to the first term and other parts of j, i.e. j^i and j^2, contribute 
to the higher order. Multiplying the above equation by A'^{y) and integrating it over y, we 
obtain: 



r J d'xd''yAlix)A^iy){T[fJix)f^riy)]) + if j d'^vKivW^^y)) = (41) 

Due to Eq.(^), the above equation is actually: 

2tf'Jd'yAl{y){r\y)) = (42) 

Surely it is the same as Eq. (|39|) itself. 



The self-energy tensor Il^'^'"-'^{x,y) = (0|T[j^o"(x)j^'o(y)]|0) corresponds to one-loop dia- 
grams of the (quantum-)gluon and the ghost with two external A legs(see fig.l). In momen- 
tum space, 11'^'^'°''' (p) is given by: 

T^lI^,abf^ . -rTl^u,ab, 



Yl^..,ab^p^ = n^"' (p) + n^'' (p), 

i(^|(^ + ln(4vr)-7 + | + ln(3#^ 



T^r\p) = i9'CA^^\g^^v' - P>^) (i^f (i + ln(4vr) - 7 + f + ln(z#3^)), (43) 

^'r\v) = ig'CA5'^\g^^v' - V^vni^kih + 1^(4^) - 7 + I + ln(-^)). 
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where Hq''"' is from the gluon loop and U'^''^ from the ghost loop; Ca = Nc for SU{Nc):, 
fi'^ is a mass scale; 7 is the Euler's constant; The pole term l/e,. = 2 — D/2 arises in the 
regularization. We calculate IIq^'" with the gauge parameter a = 1, i.e. we choose the 
Feynman gauge. The term — + ln(47r) — 7 is to be cancelled by the counterterm in MS 
scheme. After the reduction of the divergent part for n'^'^'"*, we obtain: 

W^^'-'ip) = rg^^^lcA5''\9'''p' - pV) ln(^^). (44) 

[At^Y 6 — p^ — ze 

Then we take the Fourier transformation and get 11^'^'"'' (a;,?/): 



W^''\x,y) = J ^W^^^p) exp{-ip{x - y)) 

= -tg'CAj^f6''\g'"'d'' - d^^d^H^x - y). 



(45) 



where d acts on a; or y and H{x) is the Fourier transformation of ln( _ ^_j^ ) and its derivation 
is given in Appendix B. Substituting Eq.(^5|) into Eq.(p7|), we derive: 

{f^'^^x)) = -g'CAj^,^ J d'yH{x - y){g^''d' - d^d^Aliy) (46) 

where d acts on Aiy). 

It is interesting to compare Eq.(0) with the self-energy tensor of the quark loop with 
two external A legs: 

nrC = -'3'j^\nf5'^\g'^p' - V^f) In(z^) (47) 



Note that Eq. (|4^) and (^Tf ) have opposite signs. This is just the feature of the non-Abelian 
gauge field which is of the same origin as the asymptotic freedom and the confinement. It 
means that the current generated by the gluons is opposite in sign to that by the quark. In 
Eq. (PI) , we require p^ < to give a real current j(x). If p^ > 0, the current j{x) has an 
imaginary part which would be related to the production rate for the real gluon pair. 

In fact, in the derivation of 11^^'"'' (p), there is a small positive number ie in the gluon 

2 

and ghost propagator. It finally appears in In _ ^_. . Then we can get the real part of 
YiiJ-i^,ab ^p-^ which corresponds to the production rate by noting Jm In _ ^_. = — vr. Similar to 
the definition of the production rate for qq in WM, we obtain the production rate for gg: 



Rgg = Rejd^pAl{p)W'''^\p)A^^{~p) 

= Y-^CAasjd''pe{p^)[p\A-{p)-A-{-p)) - {p-A'^{p)){p-A'^{-p))] 



(48) 



Comparing this rate with that for producing qq in the classical field: 



R,-, = -i^n;a, |dV(/)[/(^"(p) ■ A^{-p)) - {p . A\p)){p ■ A\-p))] (49) 

Note that Rqq has the same form and sign as Re+e-, the rate for the electron-positron pair 
production in QED, while Rgg has opposite sign. We remember that if A^ —>■ gg and 
Ad -^ qq are real physical processes, Rgg and Rgg can be written as Rgg = \M{Aci -^ gg)\^ 
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and Rqq = \M{Aci — > qq)]"^ respectively, therefore we have Rgg, Rq-q > 0, but from the above 
result Eg . (^ , ^91) they have opposite signs. Hence, the leading order term in g does not give 
the correct result for the gg pair production. We note that the leading order (in g) gluon 
pair production term is not gauge invariant with respect to the gauge transformation given 
by Eq.(^. In order to obtain correct gauge invariant result for gluon pair production we 
have to add higher order terms in g, which contains three and four A's, into this result. For 
details about this topic, see |]TB . 



Now we try to understand the physical implication of the current induced by the gluon. 
In the lowest order, Eq.(^) can be written as: 

D^M{^)]F''''[M^)] = -^'^^(i^y /^'^^(^ - y^^a'^d' - d^dnA^iy) (50) 

In the lowest order, this equation of motion for A corresponds to the effective action: 

Seff[A] = -Id'x\F/^Mi^)]F^''^^[A{x)] 

+9'CAj^f I d'xd''yAl{x)H{x - y){g^^^d^ - 9^9-)A^(y) ^^^> 

We can write 5*6// [A] in momentum space as: 

Seff[A] = -yd'pA-i-p)A-{p)ig^^-p^-p^^p'^) 

■[l+9'CAj^fH^J] (52) 

+ square and cubic terms of A{p) 

One can define a new amplitude A'{p) with the radiation correction caused by the gluon in 
the lowest order: A'(p) = [1 + K(p)Y^'^A(p) where k(p) = (7^T^;^yCAln( _ ^_^^ ). Therefore, 
the correction effect brought by the gluon fluctuation on the classical background field is 
similar to the medium screening effect. This screening effect becomes larger as \p'^\ decreases 
and lower as |p^| increases. This is just the feature of the effective coupling constant. In 
fact, since there is a relation between the renormalization constant Za for A and Zg for 

— 1/2 

g: Zg = Zj^ , gA should be unchanged after redefinition of the field and the coupling 
constant. Now we redefine the field as A'{p) = [1 + k(p)Y''^ A(p) , accordingly we can redefine 
the coupling constant as g^jfip) = [1 + ^(p)]"^fl'^(p)- The above argument is made for the 
gluon situation. If we include the quark sector into the classical action, we should add a term 
to k{p): k{p) = (7^74:^ [yC^ — |%] ln( _ ^_^J . It shows the screening effect brought by the 
quark is opposite to that by the gluon. We see that the effective coupling constant glffij)) 
has a running property and shows the feature of asymptotic freedom. Expressed by the new 
field and the effective coupling constant (^e//, we can write the original field equation (|50D 
as: D']^[A' , g^f f\F^^'^[A' , gef f] = 0. Therefore, the current induced by the quantum gluon 
is just the 'displacement' current in the polarized vacuum, and its effect is equivalent to 
redefinition of the field and the coupling constant. 

VI. SUMMARY 

In summary, for a system which contain a classical chromofield and the gluon simulta- 
neously, we derive the equation of motion for the classical chromofield and for the gluon. 
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Furthermore, we prove that the equation for the classical field is the same as for the gluon 
in a functional path integral sense. The proof is not straightforward, it needs a symmetry 
property between the gauge fixing and ghost sectors of the full action. The inhomogeneous 
field equation contains a current term which is the expectation value of a composite opera- 
tor including linear, square and cubic terms of the quantum field Q. We call the current as 
induced by the gluon. Note that this current essentially is classical quantity. We also derive 
the identity for this current. These identities represent the constraint on this current from 
the gauge invariance. We calculate the current at the leading order. In momentum space, 
it is just the contraction of the classical field A{p) with the self energy tensor n(j9) with two 
A's. The current generated by the gluon is opposite in sign to that generated by the quark. 
This is just the feature of the non-Abelian gauge field theory which is asymptotic free. We 
require p^ < to give a real current i{x). If p^ > 0, the current j{x) has an imaginary 
part which would give the production rate for the real gluon pair. Physically, the current 
induced by the quantum gluon is just the 'displacement' current in the polarized vacuum, 
and its effect is equivalent to redefinition of the field and the coupling constant. 

We emphasize that all the analysis done here is with respect to QCD vacuum where no 
medium effect is included yet. Normally one works in closed-time path integral formalism 
in order to derive a current in a medium. Then the current can be related to the non- 
equilibrium distribution function of the gluon. We plan to address this issue in the future 
where we wish to combine closed-time path integral formalism with background field method 
of QCD. 
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FIGURES 



A A / \ A 



FIG. 1. The one-loop diagrams for the self-energy with two external A legs. One diagram is 
the gluon loop(screw line) and the other for the ghost loop(dashed line). 



Appendix A. 

In case II, since all sources are functionals of two independent variables A and Y, the 
Ward identity (^) changes to: 

where 5W/5A^^ and 5W/SY^ are given by: 

|| = Di''[A]F'^'^'''[A] - (j^'^-) + YJ^{^) 

fi jj, fj, /J, 

Since all the field averages are set to zero during our derivation, j~ is the same as Eq.(p7|) 

and ^ = {J^''''). So we also get Eq.(|5g). 
Appendix B. 
To derive the Fourier transform of ln[— p^ — ie], we set: 

^('^^^)^l^^'"'''^''^-Po + ^"\P\"-'']- (55) 

To calculate I{x,^), we first evaluate the following function J(x,^): 

2^J{x,^)- ^^ -7(2^)4^ _p2 + ^2\j^2 _ ,, (56) 

First integrate over po- 



J(x,0 



(27r)3 J (27r) (po - ^ + *e)(po + m " ^e) 



, ^3„ „-ii\p\x°+ip-x I. jZ ii\p\x°+ip-x 
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Then integrate angular variables: 



2CJ{x,O = 0{x^ 



Q r°° p'^dp e 



-i^px^ 



+e(-x'' 



(27r)2 |f| 
p'^dp e 



r^ip\x\ _ ^-ip\x\ 



i^px" 



(27r)2 \x\ 



,^ip\x\ _ ^-ip\x\ 



Next, integrate with ^: 



[X 



0\ roo rPrlrt ^-i^P^" 



p dp e 



x\ Jo (27r)2 —ipx^ 



,^ip\x\ _ g-ip|xh 



I 0{-^') f^ P^e^^^.,|.T _ ^_.,|.^^ 



|x| JO (27r)2 ipa;0 

_e{xj_d_ r°° dp i^p^o, ip\s\ -ip\x\^ 
\x\{x^yd^Jo (27r)2 ^ ' 

6{-x^) d f^ dp 



\x\{x^ydUo (27r) 



_^iipx° /ip\x\ _ g-^plxj^ 



— « r f/ Uc 



4^(a;")2|f 



27r2^|f|(x")2 [|fP-e(a:")' + ie 



-a; 



4^(x°)2|f 



Uc uc 



0^2 



f|2-^2(a:0)2-ie] 



The final result becomes: 



lix,0 



Taking ^=1, we get: 



n^ 



(Xo) 



Oi-xo] 



f|2-e(a;")2 + ieJ2 [|f p _ ^2(^,0)2 _ ^^ 



d'^p 



^e-^P--\n[-pl + \p\'-ze] 



TT 



la^oj 



\x\ 



a;0)2 + it] 



+ 



-a^oj 



Ixp — fx^l^ — -^gl2 



Then H{x) is given by: 



H(x) 



d'^p 



-ip-x 



ln[ 



/i^ 



-p^ — ?e 



5(a;) luyU^ 



vr^ 



la:oj 

„0A2 



xr — (x^)^ + ?er \\xr — [x^Y — ie 



0\2 



(58) 



(59) 



(60) 



(61) 



(62) 
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